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Abstract. We consider the stochastic propagation of high-energy protons and nuclei in the cosmological microwave and in- 
frared backgrounds, using revised photonuclear cross-sections and following primary and secondary nuclei in the full 2D 
nuclear chart. We confirm earlier results showing that the high-energy data can be fit with a pure proton extragalactic cosmic 
ray (EGCR) component if the source spectrum is oc E~ 16 . In this case the ankle in the CR spectrum may be interpreted as 
a pair-production dip associated with the propagation. We show that when heavier nuclei are included in the source with a 
composition similar to that of Galactic cosmic-rays (GCRs), the pair-production dip is not present unless the proton fraction 
is higher than 85%. In the mixed composition case, the ankle recovers the past interpretation as the transition from GCRs to 
EGCRs and the highest energy data can be explained by a harder source spectrum oc E~~ 12 - E~ 23 , reminiscent of relativistic 
shock acceleration predictions, and in good agreement with the GCR data at low-energy and holistic scenarios. 
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1. Introduction 

One of the keys to understanding the origin of cosmic-rays 
(CR) is the spectral shape of the transition from Galactic 
cosmic-rays (GCR) to extragalactic cosmic rays (EGCR). The 
Galactic origin of low-energy CRs is generally accepted, while 
the highest energy CRs (above ~ 10 19 eV) are no longer con- 
fined by Galactic magnetic fields and most probably originate 
from other galaxies (see, however, Dar and Plaga, 1999; Blasi 
et al., 2000; Plaga, 2002). Thus, a transition between the two 
components has to occur in some energy range. The most natu- 
ral location for this transition is around ~ 3 10 18 eV at a feature 
in the CR spectrum known as the ankle. This is the only en- 
ergy range where the spectrum gets harder (i.e., its logarithmic 
slope gets smaller), offering a simple transition scheme to a 
harder CR component which is subdominant at lower energies. 

A different conclusion has recently been proposed on the 
basis of composition results from the HiRes Collaboration, ten- 
tatively showing a transition from heavy to light primary nu- 
clei at an energy around 5 10 17 eV (Abbasi et al. 2005), which 
can be identified with a second knee feature in the spectrum. 
Such composition measurements rely on statistical determina- 
tions of the elongation rate of extensive air-showers and depend 
on shower development simulations which are still quite uncer- 
tain. One should thus remain cautious about such a result, but 
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taken at face value it suggests a transition from a heavy GCR 
component to a light EGCR one at an energy below the ankle. 
A phenomenological class of models where the GCR/EGCR 
transition is around the second knee can be designed such that 
the extragalactic component accounts for CRs down to ener- 
gies ^ 10 18 eV (Berezinsky et al., 2004). In these models, the 
ankle is interpreted as an e + -e~ pair-production dip resulting 
from proton propagation over large distances in the cosmic mi- 
crowave background (CMB). This effect is the pair production 
analog of the GZK flux suppression due to photo-pion pro- 
duction at higher proton energies (Greisen, 1966; Zatsepin and 
Kuzmin, 1966). As in any model where the spectrum changes 
from a harder to a softer component, this class of models re- 
quires some unavoidable "fine tuning". Part of the fine tuning 
in the region of a GCR/EGCR transition can be explained if the 
extragalactic magnetic field has a high enough (but reasonable) 
intensity to prevent lower energy CRs from reaching us from 
distant galaxies (Lemoine, 2004; Aloiso and Berezinsky, 2004; 
see also Parizot, 2004, for a comparison of the critical magnetic 
field with an equipartition value). While the composition mea- 
surements and the matching of GCRs to EGCRs at the second 
knee remain unclear, the main point to note is that the transition 
does not have to be at the ankle. 

Phenomenological approaches are useful since none of the 
GCR and EGCR components are well understood yet, and their 
origin remains unknown (see, e.g., Parizot, 2005). Concerning 
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the EGCR component, although major uncertainties remain in 
the ultra-high-energy (UHE) range, where the expected GZK 
suppression of the flux may not be found in the data (e.g., the 
AGASA spectrum in Teshima et al., 2003), propagation studies 
find that the measured spectrum is best fitted by assuming a 
power-law source spectrum in E~@ with f3 - 2.6 (De Marco et 
al., 2003), as also found by the second-knee-transition models. 
These studies, however, assume that EGCRs are only protons. 

In this Letter, we investigate how these conclusions are 
modified when nuclei heavier than protons are also present 
in the source. We apply the revised scheme of photo-nuclear 
interactions described in Khan et al. (2005). Interestingly, we 
find that the two main aspects of EGCR phenomenology are 
jointly affected, namely the best fit slope of the source spec- 
trum and the interpretation of the ankle. The interpretation of 
the ankle as a pair production dip is not viable in the mixed 
composition case, unless the fraction of protons is above 85% 
and spectral indices are fine tuned to ad hoc values. When rea- 
sonable assumptions for the fraction of nuclei in the EGCR are 
introduced, the traditional interpretation of a GCR/EGCR tran- 
sition at the ankle is recovered and the best fit spectral index 
for EGCR is ~ 2.2 - 2.3, as expected from relativistic shock 
acceleration. Finally, the best fit EGCR spectrum occurs when 
both the GCR and EGCR component have the same spectral 
index, allowing for holistic approaches to the CR spectrum. 

2. Phenomenological source model 

Although extragalactic magnetic fields could modify the 
UHECR spectrum (e.g., Medina Tanco, 2001; Isola et al. 2002; 
Deligny et al., 2004), we assume here that their intensity is 
low enough to be negligible and focus on a comparison be- 
tween pure proton and mixed-composition models. The choice 
of a source composition is arbitrary in the absence of a source 
model, but we consider a "generic composition" assuming that 
EGCRs have the same relative source abundances as the best 
known, low-energy GCRs, derived from Ulysses and HEAO 3 
data (Duvernois and Thayer, 1996). We convert the differential 
abundance ratios, x,, at a given energy per nucleon (E/A) into 
ratios at a given energy, £ = XiAf^ 1 , where a is the spectral 
index, so that the GCR source spectrum is: Ni(E) oc ^E^ a . 

Another ingredient of our source model is the spectral in- 
dex at high energy, /3, which may be different from a, raising 
the difficulty of whether we should use x, A?~ 1 E~P or x, Af 1 E~& 
at the source. Since the CR composition depends essentially 
on the injection mechanism and the latter is uncertain for both 
EGCRs and GCRs, there is no reason, in principle, to expect 
that both components have the same differential composition. 
On the other hand, if EGCRs are accelerated out of the inter- 
stellar medium by a similar type of plasma process as GCRs, 
it is not unreasonable to assume a similar composition as well. 
Since our main goal is to investigate the difference between 
pure proton and mixed CR compositions, we study the general 
case of XjA"~ l E~P and a = ft as a special case. Although we ex- 
plore j3 between 1 and 2.7, physically motived spectral indices 
are typically (3 ~2.2 - 2.3 for the source spectrum. 

Finally, we choose a rigidity dependent maximum energy at 
the source for the various nuclear species E m . dx (^X). We assume 



that energy losses and photo-fragmentation inside the source 
can be neglected, so that all nuclei with the same gyroradius 
behave the same way. This implies E maK (^X) = Z X £ max (}H). 

3. Propagation model 

We compute the propagated spectrum of EGCRs for a uni- 
form source distribution with negligible magnetic fields, using 
a Monte-Carlo technique. For protons, we take into account the 
energy losses due to the photo-production of pions and e + -e~ 
pairs, caused by the interaction with CMB photons. Photo-pair 
production is treated as a continuous process, as allowed by the 
short interaction length and small inelasticity, while pion pro- 
duction is treated stochastically: we first draw randomly the en- 
ergy of the interacting photon in the proton rest frame, from the 
relevant boosted black-body distribution weighted by the inter- 
action cross-section, and then calculate the inelasticity from the 
reaction kinematics after choosing randomly the interaction an- 
gle in the center-of-mass frame. We also take into account the 
CMB temperature change as a function of time, as well as the 
energy losses due to the expansion of the universe. 

The propagation of nuclei is also followed stochastically. In 
addition to e + -e~ pair production, nuclei are subject to photo- 
erosion, i.e., they lose nucleons through photo-nuclear interac- 
tions with the CMB and infrared background (IRB). We use 
here for the first time a 2D scheme to follow the nuclei in the 
(A,Z) space, thanks to both revised and new cross-sections as 
described in Khan et al. (2005). This includes the following 
processes: i) the giant dipolar resonance (GDR), with a loss of 
one or more nucleons, as well as a-particles (this process oc- 
curs with photons above a threshold of ~ 8 MeV in the nucleus 
rest frame), ii) the quasi-deuteron (QD) process, where a vir- 
tual pion interacts with a nucleon pair within the nucleus, lead- 
ing to the ejection of the pair and possibly additional protons 
or neutrons (this occurs with photons of typically 20 MeV), 
iii) baryonic resonances (BR), where a real pion is produced, 
ejecting a nucleon and possibly interacting further with a nu- 
cleon pair, eventually leading to the loss of an average of six 
nucleons for an Fe nucleus (this occurs above ~ 150 MeV), 
and iv) the photo-fragmentation (PF), occurring at very high 
energy (~ 1 GeV) and breaking the nucleus into many frag- 
ments of much lower mass and energy. 

The new GDR cross-sections are discussed in detail in 
Khan et al. (2005). For the higher-energy processes, we used 
the parameterisation of Rachen (1996). The QD cross-section 
is taken as o-Q D (e,A) = (0.55 mbarn) A 5/4 (e/eQ D - l) 3/2 e/e Q Q, 
where e is the photon energy in the nucleus rest frame, A is 
the nucleon number and eQ D ^ 33.3 A~ 1/6 MeV is the in- 
teraction threshold. The number of emitted nucleons is ob- 
tained from the experimental branching ratios of Puget et al. 
(1976). The BR cross-section is modeled as the sum of four 
Lorentzian representing the four main resonances, cr BR (e,A) = 
A Yn=\ CjClC^ ; Ei, T,), with the following parameters (in GeV, 
GeV and //barn, respectively): (£,-, F,-, cr,) = (0.34,0.17,351), 
(0.75, 0.50, 159), (1.0, 0.60, 21) and (1.5, 0.80, 26). Finally, the 
PF cross-section is modeled as cr PF (e, A) = (1 jifbarn) A 91 (l - 
2.48 exp(-e/0.8GeV))(69.8s 081 + 64.3 s' 0453 ), where s = 
0.88 + e/(532MeV). We also follow Rachen (1996) for the 
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Fig. 1. Left: mean-free-path of 56 Fe nuclei against the various photo-nuclear processes in the CMB and IR background (see text), as a function of 
T. Right: energy loss length (including pair production) for various nuclei (labeled by their original A) as a function of E, with their dispersion. 



pair production cross-section, including corrections to the for- 
mulae of Blumenthal (1970). The attenuation time of a nu- 
cleus with Lorentz factor F is related to that of the proton, t p , 
as: t(Z, A, y) = (A/Z 2 ) r p x <3> _1 (Za), where a is the fine struc- 
ture constant and d>(x) = 1 - 0.29x 2 + 0.25x 4 - 0.25x 6 . 

The propagation of all the nuclei (distributed in energy ac- 
cording to the source spectrum) is simulated with a Monte- 
Carlo code taking the above processes into account. We keep 
track of all secondary nuclei emitted and propagate them in the 
same way. It is assumed that each ejected fragment (or nucleon) 
takes away the same energy per nucleon as the parent nucleus, 
so that it has the same Lorentz factor. In addition to the CMB, 
we include the IRB as described in Malkan and Stecker (1998): 
n{E y ) =1.1 \Q- A E- 15 for 0.02 eV < E y < 0.8 eV. The result- 
ing mean free path for a nucleus with Lorentz factor F is given 
by: 

AiR 1 = deA T e'cr(e')de', (D 

11 Jf„,/2r e Je, h 

where e t h is the energy threshold of the process under consid- 
eration and e max = 0.8 eV. 

The total mean free path of 56 Fe nuclei against the vari- 
ous photo-nuclear processes in the CMB and IRB is shown in 
Fig. la as a function of F. The GDR process clearly dominates 
up to F = 10 10 5 — 10 1 1 and is responsible for the main GZK 
effect for nuclei (sudden reduction of the flux). In Fig. lb, we 
show the length scale of energy losses for protons, He, O and Fe 
nuclei at the source, defined as the propagation length, xis {E), 
after which a nucleus of initial energy E has lost 25% of its en- 
ergy (including photo-pair production). The error bars indicate 
the dispersion due to the stochastic propagation, emphasizing 
the importance of a Monte-Carlo approach. As can be seen, the 
drastic reduction of the propagation horizon, due to the GDR, 
occurs at different energies for different nuclei and is much 
more pronounced for heavy nuclei than for protons. It should 
be stressed that, while the GZK effect appears stronger for nu- 
clei than for protons, nuclei photoerosion products contribute 
to the EGCR flux at lower energies and substantially modify 
the shape of the spectrum. This leads to important changes in 
the global phenomenology of EGCRs, as we now discuss. 



4. Results and discussion 

In Fig. 2a, we show the resulting propagated spectra ob- 
tained with the prescription described above where ^max^H) = 
3 10 20 eV. The error bars represent the l-cr fluctuations of the 
flux in each energy bin, as expected for a data set with the statis- 
tics of the AGASA experiment, i.e., 866 events above 10 19 eV 
(cf. De Marco et al., 2004). We compare our results with the 
data from HiRes (Abbusi et al. 2004) and AGASA (Takeda et 
al. 2003) experiments. In the pure proton case, we find again 
that the best fit spectrum has f3 = 2.6 and a pair production dip, 
as in previous work. However, if the composition is mixed, the 
choice of a — f3 — 2.6 (corresponding to 24% of protons at 
the source above 10 18 eV) provides a very poor fit to the data 
(cf. Fig. 2b, lower curve). Figure 2b shows how the spectrum 
around the ankle and the GZK bump evolves with the fraction 
of protons in the source composition (above 10 18 eV). For frac- 
tions lower than 75%, the pair production dip is strongly atten- 
uated and the ankle is actually reversed for fractions lower than 
~ 60%. Therefore, the interpretation of the ankle as an e + -e~ 
pair production dip requires a very large fraction of protons at 
the source (see also Berezinsky et al., 2005). In contrast, Fig. 2a 
shows two good fits to the highest energy data using mixed 
composition models with a = f3 — 2.2 or 2.3, corresponding 
to 50% or 40% protons at the source, respectively. Note that 
the observed composition at Earth is much lighter than at the 
source, because of the massive production of secondary nucle- 
ons peaking around 5 10 18 eV and the successive cut-offs of 
nuclei with increasing mass (Allard 2004). 

As can be seen in Fig. 2a, the pure proton model with 
(3-2.6 provides a good fit to the data down to ~ 10 18 eV 
(in a scenario with uniform source distribution and no mag- 
netic field). In such a model, the transition from GCR to EGCR 
should thus occur at lower energies around the second knee 
(Berezinsky et al., 2004). However, no acceleration process is 
known to yield a spectrum as steep as (3 — 2.6. Instead, both 
analytical and numerical studies of particle acceleration in rela- 
tivistic shocks give injection spectra with f3 2.2-2.3 (Bednarz 
and Ostrowski, 1998; Kirk et al., 2000; Lemoine and Pelletier, 
2003). As shown in Fig. 2a, such spectral indices can provide 
equally good fits to the high energy data, assuming a realistic 
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Fig. 2. Left: Propagated spectra for sources with protons only and ji = 2.6, and with mixed compositions and /? = a = 2.2 and 2.3. Data points 
from the main experiments are indicated, with an arbitrary downward shift of AGASA's energy scale by 20%. Right: evolution of the spectrum 
with the fraction of protons at the source, for an E~ 2 6 source spectrum. 



source model with a mixed composition a [3. In such a sce- 
nario, the transition from GCRs to EGCRs occurs at the ankle, 
which thus keeps its "standard" interpretation. 

In conclusion, our calculations show that the phenomeno- 
logical interpretation of the CR data at high energy strongly 
depends on the source composition. We confirm that pure pro- 
ton sources require a steep spectrum in ~ E~ 2 6 and are able 
to account for the ankle down to 10 18 eV, interpreted as a pair 
production dip of the EGCR component. On the other hand, if 
the EGCRs are accelerated out of the interstellar medium with 
a source composition roughly similar to that of the GCRs, we 
find that a harder source spectrum ~ E~ 2 ' 2 gives an equally 
good fit to the data, allowing an interpretation of the ankle as 
the transition from GCRs to EGCRs, with the GCR component 
extending up to ^ 10 18 eV. 

These results are interesting for several reasons. First, a 
source spectrum with /3 2.2-2.3 appears quite plausible from 
the theoretical point of view (relativistic shock acceleration). 
Second, a source spectrum ~ E~ 2 3 has also been considered as 
the best fit to the low-energy data for the GCR component (e.g. 
Strong and Moskalenko, 2001; Ptuskin, 1997). It is thus partic- 
ularly interesting to note that assuming a similar composition 
for both the GCR and EGCR components is precisely what also 
makes a similar source spectrum possible. In addition, it has 
been shown in Parizot (2005) that a source power-law index of 
~ 2.3 is a necessary condition for holistic models, in which the 
same sources produce the CRs at all energies. The results pre- 
sented here show that a self-consistent model can be built for 
CRs with a similar composition and spectrum at all energies 
provided the GCR/EGCR transition is at the ankle. Finally, the 
inclusion of nuclei with a standard CR composition at high- 
energy can avoid the fine-tuning necessary to match the GCR 
and EGCR components in the second-knee transition models. 
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